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The nonlinear propagation of electromagnetic (EM) electron-cyclotron waves (whistlers) along 
an external magnetic field, and their modulation by electrostatic small but finite amplitude ion- 
acoustic density perturbations are investigated in a uniform quantum plasma with intrinsic spin of 
electrons. The effects of the quantum force associated with the Bohm potential and the combined 
effects of the classical as well as the spin-induced ponderomotive forces (CPF and SPF respectively) 
are taken into consideration. The latter modify the local plasma density in a self-consistent manner. 
The coupled modes of wave propagation is shown to be governed by a modified set of nonlinear 
Schrodinger-Boussinesq-like equations which admit exact solutions in form of stationary localized 
envelopes. Numerical simulation reveals the existence of large-scale density fluctuations that are 
self-consistently created by the localized whistlers in a strongly magnetized high density plasma. 
The conditions for the modulational instability (MI) and the value of its growth rate are obtained. 
Possible applications of our results, e.g., in strongly magnetized dense plasmas and in the next 
generation laser-solid density plasma interaction experiments are discussed. 

PACS numbers: 52.27.Aj, 52.35.Hr, 52.35.Mw 



I. INTRODUCTION 

Having been discovered more than a century ago 
whistler waves become one of the most important waves 
in plasmas. Such waves (also known as helicons in 
solid state plasmas) are low-frequency (If) (in compari- 
son with the electron-cyclotron frequency, ui c ) right-hand 
circularly polarized (RCP) electromagnetic (EM) waves 
guided almost along the external magnetic field in dense 
plasmas. Because of the increase of their group velocity 
with the frequency, u < uj c /2 (see, e.g. Ref. .2]), the If 
waves arrive later giving rise a whistling down-effect ob- 
served at ground level. Stenzel in his classic paper Q 
demonstrated experimentally the creation of magnetic 
field-aligned density perturbations excited by the pon- 
deromotive force exerted by the EM whistlers. 

Whistler waves are important not only in space plas- 
mas due to wave-particle interactions, but also in labora- 
tory plasmas as helicons for efficient plasma production 
as well as in dense astrophysical environments (3-floj . On 
the other hand, large amplitude whistlers propagating in 
a magnetized plasma can initiate a great variety of non- 
linear effects, e.g., three- wave interactions, parametric 
instabilities [10|, modulational instability and the subse- 



quent soliton formation [J-|6( . The latter which, in turn, 
causes local electron density enhancement or depletion 
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in plasmas, are considered as a basis for understanding 
laser energy deposition in pellets 11], pulsar radiation in- 
teraction with the ambient magnetosphere [l2j ]. whistler 
wave propagation in solar winds pj$[ etc. Recent labo- 
ratory experiment [3] and observations from the Freja 
satellite [15[ show the clear evidence for the formation of 
whistler envelope solitons accompanied by plasma den- 
sity cavities. Moreover, electrons in Van Allen radia- 
tion belts can be accelerated to MeV energies within 
a short period by large amplitude whistlers [r|. The 
latter have recently been observed by the Cluster space- 
craft 



17], the STEREOS [16j and the THEMIS \tt 



Furthermore, laboratory experiments and theoreti- 
cal confirmation [2(J have demonstrated the existence of 
propagating whistler spheromaks with fields exceeding 
the ambient magnetic field. Whistlers also contribute to 
fast magnetic reconnection and plasma dynamics in two- 
beam laser-solid density plasma interaction experiments 
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Recently, there has been a notably growing interest in 
investigating various quantum plasma effects in view of 
some experimental progresses in nanoscale plasmas [23| , 
ultracold plasmas 24[, spintronics [25| and plasmonics 
26]. On the other hand, superdense quantum plasmas 
are omnipresent in compact astrophysical objects, e.g., 
the interior of massive white dwarfs, interior of Jupitors, 
magnetars etc. [27I - I29I ] . as well as in the next genera- 
tion intense laser-solid density plasma interaction exper- 
iments 13- 32j. In dense plasmas, degenerate electrons 
follow Fermi-Dirac pressure law, and there are typically 
quantum force associated with the Bohm de Broglie po- 
tential, which produce wave dispersion at nanoscales [33l — 
[35j ] - Furthermore, the effects of the electron spin man- 
ifests itself in terms of a magnetic dipole force, as well 
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spin precession, which can be exploited by transforming 



More 



the Pauli equation to fluid-like variables 
elaborate kinetic models has also been developed |38l.l39|. 
Hence the dynamics of electrons in Fermi degenerate 
plasmas will be affected not only by the Lorentz force, 
but also by the effects of quantum statistical pressure, 
the Bohm force as well as the effects due to intrinsic spin 
of electrons. We ought to mention that in a dense mag- 
netized plasma there also exist spin waves, which can 
be excited by intense neutrino fluxes. Thus, nonlinear 
theories of EM waves, in particular whistlers in magne- 
tized dense plasmas need to be developed in its own right 
accounting for all these quantum effects. Recently, the 
theory of the ponderomotive force in plasmas has been 
extended to account for the contribution from the in- 
trinsic spin of electrons [loj]. It has been demonstrated 
that an EM pulse can induce a spin-polarized plasma by 
this spin-ponderomotive force (SPF). Such force could 
also play an important role in the propagation of If EM 
waves, e.g., whistlers, Alfven waves. 

Our objective here is to present a theoretical study 
of modulated whistler wave packets interacting nonlin- 
early with background If density perturbations that are 
reinforced by the classical ponderomotive force (CPF) [f| 
as well as the SPF 40]. The role of the ion motion as 
well as the dispersive effects due to charge separation and 
the electron tunneling are also taken into account. We 
will include the field aligned velocity perturbation (free 
electron streaming) associated with the If motion, and in 
addition, generalize the related classical results that ex- 
ist in the literature (see, e.g., Refs. 0,0]). The obtained 
results could be useful for understanding the propagation 
of localized EM whistlers which may emanate in the in- 
terior of magnetized white dwarfs, magnetars as well as 
in the next generation intense laser-solid density plasma 
experiments. 



II. NONLINEAR EVOLUTION EQUATIONS 

Let us consider the propagation of nonlinearly coupled 
EM whistlers and ion-acoustic (IA) density perturbations 
along a constant magnetic field B = Bqz in a quantum 
electron-ion plasma where any equilibrium drift velocity 
is zero. In the modulational representation, the high- 
frequency (hf) EM wave field for the RCP whistlers is 
given by E = (x — iy) E(z, t) exp(ikz — iujt)+c.c, where 
E(z,t) is the slowly varying (both in space and time) 
envelope of the whistler wave electric field and c.c. stands 
for the complex conjugate. Also, u) (k) represents the 
whistler wave frequency (number). The basic equations 



for the evolution of nonlinear whistlers then read 36 , 4(| 
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where n e , m e , v e denote the number density, mass and 
velocity of electrons respectively, B is the magnetic field 
and P e is the electron thermal pressure. Also, S is the 
spin angular momentum with its absolute value |S| = 
| So | = ft/2; M = - (ff/2) Mb, where g 2.0023193 is the 
electron g-factor and /is = eh/2m e is the Bohr magne- 
ton. The equations ([I])-© are then closed by the follow- 
ing Maxwell equations with V.B = 0. 
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The equations CO)-© represent the nonrelativistic evolu- 
tion of spin— 1/2 electrons, and are applicable even when 
different states with spin-up and spin-down (relative to 
the magnetic field) can be well represented by a macro- 
scopic average. This may, however, occur in the regimes 
of very strong magnetic fields (or a very low temperature 
regimes) , where generally the electrons occupy the lowest 
energy spin states. On the other hand, for a time-scale 
longer than the spin-flip frequency, the macroscopic spin 
state is well-described by the thermodynamic equilibrium 
spin configuration, and in this case the above model can 
still be applied. However, such case in which the macro- 
scopic spin state will be attenuated by a factor decreasing 
the effective value of |S| below H/2, will not be consid- 
ered further in the present work . As a consequence, our 
studies will be focused on the regime of strong magnetic 
fields and high density plasmas. 

Taking the curl of Eq. © and using Eqs. ©-(0) 
we readily obtain the following evolution equation for 
whistlers. 
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In the linear theory, the whistler frequency uj and the 
wave number k are related by the following linear disper- 
sion relation in the nonrelativistic limit (see for details, 
Ref. 1). 
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where = ck/u is the refractive index, uj^ = 

g 2 \Sq \ /Am e X 2 is the frequency due to the plasma magne- 
tization current and A e = c/uj pe is the electron skin depth 
with w pe(i) 



y/n e 2 /e m e ^ 
plasma frequency. Also, uj c 
cyclotron frequency and w g 
spin-precession frequency. 



denoting the electron (ion) 
= eBa/m e is the electron- 
= (g/2)uj c is the electron 



The nonlinear dynamics of whistler wave envelopes un- 
der the modulation of electron density perturbations as- 
sociated with the If IA fluctuations and of the nonlinear 
frequency-shift caused by the magnetic field aligned free 
streaming of electrons with flow speed w ez , can be de- 
scribed by the following nonlinear Schrodinger (NLS)-like 
equation which is obtained from the EM wave equation 
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where E = E x — iE y , and the group speed, v g = dui/dk 
[see Eq. (11) in Ref. 0] and the group dispersion, v' g = 
d 2 uj I dk 2 of whistlers are given by 
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The nonlinear frequency shift A is given by 
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where A = 2c 2 k/w 2 e + g 2 hk/4m e (uj — uj g ) and N = 
n e /no is the relative perturbed density. By disregard- 
ing the spin contribution one can recover the previous 
results 0, HI • Note that the term oc v ez , representing the 
Dopplcr shift due to the plasma streaming along the ex- 
ternal magnetic field, is no longer negligible, but may be 
comparable to the other nonlinear terms, and can thus 
change the sign of the nonlinearity as well. More pre- 
cisely, both v' g and A will change their sign depending 
on the frequency range to be considered as well as the 
contribution from the spin correction terms. Later, we 



will see that the change of sign is important for the for- 
mation of localized wave packets at different whistler fre- 
quencies. The quantities N and v ez are related to each 
other by the electron continuity equation. 
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Note that the ponderomotive force due to the EM 
whistlers usually drives the If (compared to the whistler 
wave frequency uj) density perturbations which propagate 
along the field lines with low-phase speed (compared to 
the electron thermal speed). Thus, the If electrostatic 
modulation also satisfies the electron momentum equa- 
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where Ei is the If part of the wave electric field and Vp 
= y/ksTF /m e is the Fermi speed relevant for a high den- 
sity plasma 42|. Here Tp = h 2 (37r 2 no) 2 ^ 3 /2fcsm e and 
ks is the Boltzmann constant. The term oc h 2 is the 
quantum correction associated with the Bohm de Broglie 
potential. The ponderomotive force contributions are 
proportional to the constants T\ and T2 where 

g 2 hk 2 
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in which the first terms appear due to CPF [5| and the 
second ones (oc K) are due to the SPF [40|. The equations 
for the cold ion motion involved in the If IA perturbations 
are 
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Eliminating rij, J5j, Vi Z and disregarding the term oc 
m e /mi, we obtain from Eqs. (fT3|) . ([T5|) - ([T7| the driven 
wave equation for If perturbations of the Boussinesq-type 
as 
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where /ii = eoTi/2nom e and /12 = £akl '2 / '2nom e uj 2 , c s — 
\JkBTp /mi is the ion-acoustic speed and Xf = c s /uj p i is 
the Fermi screening length for electrostatic oscillations. 

Thus, we have a set of three coupled equations, namely 
©, (TT^I) and (TT5)) . modified from previous results by the 
SPF and quantum tunneling, which describes the non- 
linear coupling of electron whistler waves with the field 
aligned electrostatic density fluctuations. These equa- 
tions can be recast by normalizing the variables accord- 
ing to z — > z/Xp,t ~ > tuj p i,E — > E/Eq, v ez — > v ez /c s , in 
which case we obtain 
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where E Q = y/2k B T F n /e , V g = v g /c s ,Vg = v' g uj pt /c 2 s , 
\& = A/wpj, H = hw pe /2k B T F is the quantum coupling 
parameter, Ai = w 2 e Fi/o; 2 and A2 = uj 2 e kc s T2/uJ 2 - Equa- 
tions (| 1 9 1) - (f2T|) contain the main results of the present 
work. In particular, previous results 0, H| can be recov- 
ered by disregarding the spin contribution oc h as well 
as the particle dispersion oc H and considering, e.g., the 
isothermal equation of state (relevant for low or moderate 
density plasmas). 



III. STATIONARY LOCALIZED SOLUTIONS 

In this section we will investigate the properties of non- 
linear whistlers by solving numerically the Eqs. (|19p - (f2Tj) 
in the stationary frame £ = z — Mt (where M = V/c s ). 
We will consider the parameter regimes for the density 
and the magnetic field for which the nonrelativistic fluid 
model is valid and SPF is comparable to the CPF. We 
will also see that the case in which SPF dominates over 
the CPF may correspond to the strongly magnetized su- 
perdense plasmas where relativistic treatement may be 
necessary. However, before going further to such discus- 
sions let us first consider the particular case in which 
the dispersion due to charge separtion (quasineutrality) 
is negligibe. The latter can be justified even when the 
spin effects dominate, i.e., \ = fik 2 /m e uj ^ 1 2]. From 
the scaling 



k 2 \ 2 p 
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we find that the quasineutrality limit A; 2 A|. <C 1 holds 
in nonrelativistic (Vp <§; c 2 ) plasmas as long as tir = 
ck/u > 1 and uj p i < u < uj pe . However, we will see 
that in a specific parameter regime, such restrictions can 
be valid for very If {uj <C lu c ) whistler modes. In this 
case, B < B Q = 4.4 x 10 9 T and n > 10 32 m" 3 
with Tp > Tb — huj c /kB- Moreover, when ^ » 1, the 
contribution from the term oc A2 can be smaller than that 
oc Ai, since IA2/A1I ~ (kc s /u) (m e uj/hk 2 ) <§; 1. Thus, in 
the quasineutral regime, we obtain from Eqs. (|19l) - ([2"Tj) 
the following NLS equation. 
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together with 



N = A\E\ 
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Then we can write Eq. (fTTj) as A = — A \E\ 2 , where A is 
defined as 
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where A = (A x + A 2 M) / (M 2 - l) . Physically, the elec- 
trons experience a longitudinal force exerted by the front 
of the whistler pulse, and thereby gain a net energy. The 
electrons gain energy during the rising front of the pulse, 
but then slows down by the backward ponderomotivc- 
like force. Moreover, electrons can approach the group 
velocity of the whistler when it reaches the pulse peak 
at the center. From Eq. ([23)1 . we find that this can be 
possible for M 2 <C 1, which may happen for a whistler 
frequency satisfying lo c /2 < u> < uj c and for high density 
(~ 10 36 m -3 ) and strongly magnetized (B ~ 10 8 T) plas- 
mas. In this case, the Fermi speed may exceed the group 
speed (~ c). On the other hand, corresponding to the pa- 
rameters as in Fig. 1 below, M 3> 1 and A1+A2-M ~ 0.01, 
so that N ~ \E\ 2 x 10~ 7 , v ez ~ \E\ 2 x 10~ 5 and A - 10~ 7 . 
Again, note that slow electrons can freely move along 
the direction of the external magnetic field. The finite 
velocity perturbations would then induce an additional 
density change in order to maintain the conservation of 
particles (equation of continuity) under localized distur- 
bances. Consequently, the total density variation in the 
frequency-shift becomes 5N, where 
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Clearly, 5 changes sign whenever the third term oc h in 
Eq. (|26l) dominates over the other two terms. Now, for 
If propagation of whistlers, Ai < and as in the previous 
section, A 2 Af (> 0) is smaller compared to Ai when the 
spin contribution dominates. Thus, in the quasineutral 
If regime, the density and velocity perturbations are pos- 
itive and negative according as the whistler wave propa- 
gation is subsonic or supersonic [see Eq. ([24)l ]. 

Furthermore, localized bright (dark) envelope solutions 
of Eq. (125)) exist through the modulational instability 
(stability) when Av' g > 0(< 0). For If waves (w < u c ), 
when hk 2 /m e tu > 1, Ai + A 2 M < 0, v' g < and A ^ 
according as M 1. Hence, a possible final state of the 
MI could be a supersonic (subsonic) bright (dark) soliton- 
like structure in a quasineutral spin quantum plasma. 
Equation ([25)) has an exact soliton solution (when A and 
Vg have the same sign) of the form 
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FIG. 1: (Color online) Whistler solitary solution of Eqs. Ip8)l 
and (|29p with associated electric field W (upper panel) and 
density perturbation iV (lower panel) for = 0.8 (solid line) 
and 0.82 (dashed line). The other parameter values are no — 
10 34 m^ 3 , B = 5 X 10 8 T, lj c /lu pc = 15.6, uj = 0.4, M{= 



Vg/Cs 



234.54, V' a 



-1.81, v a 



x 10 s m/s. 




FIG. 2: (Color online) Whistler solitary solution of Eqs. ((28)) 
and (|29p with associated electric field W (upper panel) and 
density perturbation N (lower panel) for Bo = 5.8 x 10 8 T 
(solid line) and 5.9 x 10 8 T (dashed line). The other parameter 
values are no = 10 34 m -3 , Q — 0.8, uj = 0.4. The correspond- 
ing M values are M = 234.88 (solid line) and 234.91 (dashed 
line) . 
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where E m ,^o are constants. The other particular cases, 
namely the quasistationary If density response (i.e., d t — > 
0) for which to « to pi Q and the case of unidirectional 
propagation (near sonic envelopeHn which the quasineu- 
trality is not a valid assumption [7fl will not be discussed 
here as those cases are not so relevant to the parameter 
regimes to be considered, instead we will focus on our 
main Eqs. (19])- (EH- 

Thus, we look for stationary solutions of Eqs. (fT9")l - 
(j2"Tj) in the stationary frame £ = z — Mt. Here we assume 
E to be of the form E = W (£) exp (— iQt) ,where W is a 
real function and is a real constant. Then Eqs. (|T9|) - 
(|2Tj) reduce to 

V d 2 W 

-J7t-+Wn + ANW = 0, (28) 
2 at;* 

(-M 2 + H 2 )^ + (M 2 -\)N 

( d 2 W 2 \ 
= (Ax + A 2 M) ( M 2 —- + W 2 ) , (29) 

where A = A/A. Wc numerically solve the equations 
(|28|) and (|29|) by Newton method with the boundary 
conditions N, W, d^N, d 2 W -> as |£| -> oo. We con- 
sider the density and magnetic field strength to vary as 
n ~ 10 34 - 10 36 m" 3 and B ~ 10 8 T. Figure 1 illus- 
trates the existence of double-hump localized whistler 
envelope accompanied with a density depletion for a 
set of parameters: n ~ 10 34 m -3 , B a ~ 5 x 10 8 T, 
cj = 0.4 and = 0.2. The corresponding frequen- 
cies are to pl = 1.32 x lO^s" 1 , to pe = 5.64 x lO^s" 1 
and lu c = 8.79 x 10 19 s _1 . Also, M(= v g /c s ) = 234.54, 
Ar> e (= Xfc = VFe/^pe) = 9.67 x 10~ 12 m and V Fe = 
5.46 x 10 7 m/s. Thus, the whistlers have negative group 
dispersion with VI — —1.81. From the dispersion rela- 
tion we obtain k = 1.18 x 10 11 m _1 , which corresponds to 
whistlers with a wavelength of 5.3121 x 10~ n m, and the 
group speed is v g — 2.98 x 10 8 m/s. Furthermore, the non- 
linear frequency shift is obtained as A = 0.85. The den- 
sity depletion is observed quite small due to large group 
velocity (compared to the sound speed) of the whistler 
waves. 

In another illustration (Fig. 2) with a higher magnetic 
field, we observe a dark-soliton-like structure correlated 
with a density hump. The amplitude of the solitary pulse 
decreases as the magnetic field increases. In Fig. 3 we 
have presented the solitary structures when the density 
is very high (no ~ 10 36 m~ 3 ). This basically corresponds 
to the case when Hk 2 /m e aj > 1. However, in this case 
one must note that the Fermi speed is close to or can 
even be larger than the speed of light in vacuum and so, 
nonrelativistic quantum fluid model may no longer be ap- 
propriate. The quantum parameter H has no significant 
role for the regime considered here, as can be seen that 
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FIG. 3: (Color online) Whistler solitary solution of Eqs. (f28|) 
and (|29[) with associated electric field W (upper panel) and 
density perturbation N (lower panel) for uj — 0.4 (solid line) 
and 0.38 (dashed and dotted ine). The other parameter values 
are B = 5 x 10 8 T, n = 10 36 m" 3 (for solid and dashed line) 
and 2.1 x 10 36 m" 3 (for dotted line), 9 = 0.2. The values of 
M are M — 35.38 (for solid and dashed line) and 22.23 (for 
dotted line). 

M (= V g ) mainly dominates in the term — M 2 + H 2 [Eq. 
(|29|)]. because of large group velocity (« c). In order that 
H can be comparable to M, one might have to consider 
relatively higher densities (> 10 36 m -3 ) and weakly mag- 
netized (<?C 10 8 T) plasmas. However, in this case the 
coefficient Ai + \%M (~ 10 5 ) will be much larger than 
the other coefficients, which might prevent any hope for 
localized solution. As shown in Fig. 4, one can excite a 
nondiverging whistler with a positive group dispersion in 
other regime, e.g., u = 0.189, 9 = 0.7, n ~ 7x 10 36 m~ 3 
and B - 5 x 10 8 T for which V' g = 0.262, H = 0.11, 
v g = 2.37 x 10 8 m/s, M = 45.24, V Fe = 2.25 x 10 8 m/s, 
c s = 5.24 x 10 6 m/s. This basically corresponds to oscil- 
latory pulse associated with a field-aligned density hump 
(A~10- 10 ). 



IV. GROWTH RATE OF INSTABILITY 

Nonlinear interaction of the hf pump EM whistlers 
(to, k) with If electrostatic field aligned perturbations 
(Q,K) gives rise upper and lower side bands with fre- 
quency and wave numbers respectively (to + Q, k + K) 
and (to — Q, k — K). The latter interacts with the pump 
and thus produces a If ponderomotive force which even- 
tually reinforces the If electrostatic oscillations. When all 
the perturbations are aligned along the external magnetic 
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FIG. 4: (Color online) Ducted whistler obtained as solu- 
tion of Eqs. (|28[) and (|29[) with associated electric field W 
(upper panel) and density perturbation N (lower panel) for 
oj = 0.189, Bo = 5 x 10 8 T, n = 7 x 10 35 m~ 3 , 9 = 0.7. The 
other parameters are M = 45.24, V' g = 0.262, H = 0.11, v a = 
2.37 x 10 8 m/s. 

field, the parametric interactions of EM waves can be de- 
scribed from Eqs. (|T9|) - (|21|l by the following dispersion 
relation. 

AK 2 Vg [K 2 (1 + H 2 K 2 ) - (1 + K 2 ) fl 2 ] 
= AV g E 2 (l - n 2 ) (Ai# + A 2 ft) (KTi + (CI) , (30) 

where C, = kuju>3/ \ De ( w — l) 2 and u, k have been nor- 
malized by w c and X^ e respectively. Some simplification 
can be in order. Note that under the quasineutrality as- 
sumption, the coefficient of CI 4 , CI 3 , and the term oc K 4 as 
well as the term oc Ai in the coefficient of CI 2 will not ap- 
pear. Also, for If propagation of whistlers (u> < ui c ) , A2 is 
smaller and thus being neglected. Moreover, the ratio of 
the term oc Q in the coefficient of CI (which appears due 
to the parallel electron streaming v ez ) and the constant 
term oc Ti scales as (k/K)oj (m e Cluj c /hk 2 ) and we need 
(k/K)u (m e Ciui c /hk 2 ) -C 1 for spin effects to be domi- 
nant. Thus, in this case the dispersion relation reduces 
to 

Cl 2 K 2 (1 + H 2 K 2 ) - Efa 2 , (31) 

where rj 2 = ±4VgAiFi/w 2 AV^ in which u>,k etc. are 
being normalized. Clearly, MI sets in for modulation 
wave numbers satisfying K\J\ + H 2 K 2 < Eorj, or K < 
K c w EqT] for highly dense medium and small K. The 
growth rate of instability (Cl = ij) is then given by 



7 a ^E 2 if -K 2 {l + H 2 K 2 ). (32) 

Hence, in the long-wavelength limit {K — > 0) maximum 
growth rate of instability can be achieved, and is roughly 
proportional to the pump wave electric field E$ and rj. For 
parameters as in Fig. 1, we obtain 7 w 2.77. It basically 
restricts the characteristic length-scale to a certain value 
for the formation of envelope solitons through MI. 

V. DISCUSSION AND CONCLUSION 

In the present investigation focusing on whistler waves 
we point out that the spin contribution is substantial 
when hk 2 /m e u! ^> 1, i.e., when hiu c /m e c 2 > 1 and 
bj 2 <C c 2 k 2 . This corresponds to the case in which the 
magnetic field strength, B > Bq and the particle den- 
sity is very high, i.e., no > 10 36 m~ 3 for which the mag- 
netic field is non-quantizing and does not affect the ther- 
modynamic properties of electrons. However, in such 
regimes, the Fermi velocity may approach or exceed the 
whistler group velocity (close to c in the present study), 
and so the nonrelativistic quantum fluid model may no 
longer be appropriate to consider. In the present work, 
we have considered Bq ~ 10 8 T and the density to vary 
in the range 10 34 < no < 10 36 m -3 in order that the non- 
relativistic fluid model is valid to some extent. Moreover 
that ojpe > u! c and the terms due to spin magnetiza- 
tion current together with the SPF are comparable to 
the classical counter parts. Furthermore, in this regime 
the velocity of electrons remains much smaller than the 
whistler group velocity (~ c). 

Since the whistler group speed is much higher than the 
IA speed, whistler solitons are not significantly affected 
by the particle dispersion associated with the Bohm po- 
tential as well as the Fermi-Dirac pressure, though the 
length scale of excitation is of the order of the Compton 
wavelength. However, those effects reduce the plasma 
characteristic wavelength of excitation. Such effects can 
be more significant in some other regimes when M < 1 
and/ or for possible excitation of the ion wakefields at 
nanoscales. Note that since degenerate electrons follow 
the Fermi-Dirac pressure law (where the Fermi tempera- 
ture is density dependent), the cold plasma limit can not 
be recovered from the present study unless one consid- 
eres, e.g., isothermal equation of state to be relevant for 
low or moderate density plasmas. Furthermore, H — > 
means that one approaches the higher density regimes 
and H = is the case when one simply disregards the 
quantum tunneling effect. 

The parameter regimes considered here can be achiev- 
able in the magnetized white dwarfs (~ 10 36 m~ 3 ) as 
well as in the next generation intense laser-solid density 
plasma experiments (~ 10 34 m~ 3 ), in x-ray free electron 
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lasers, and in plasmonic devices. One can, in princi- 
ple, go beyond the parameter regimes considered here 
(since there is no specific theoretical limit for the den- 
sity), however, we have to be careful about those param- 
eter values for the excitation of localized whistlers and for 
spin-ponderomotive force to have a role. The latter may 
dominantly accelerate the ions by separating the electric 
charges and building up a high electric field. However, 
plasma can sustain such high electric fields, and so it 
remains an attractive medium for particle acceleration, 
which is still a most important areas of research works in 
both laboratory and astrophysical plasmas. 

In conclusion, we have presented a new set of non- 
linear equations which governs the dynamics of mod- 
ulated whistlers interacting with the field- aligned elec- 
trostatic If density perturbations due to IA fluctuation, 
in a magnetized spin quantum plasma. Both the classi- 
cal as well as the spin-induced ponderomotive force has 
been considered to modify the local plasma density in a 
self-consistent manner. Numerical simulation of the gov- 
erning equations in the stationary frame [Eqs. (|28|) and 
(|29[) ] reveals the existence of supersonic stationary enve- 
lope solitons characterized by a single or double hump 
whistler wave electric fields that are trapped in a self- 
created density cavity. This happens for wave frequency 
satisfying u < co c /2 and when the whistler has negative 
group dispersion. When the whistler frequency is smaller 
than w c /4 and the group dispersion is positive at higher 
densities, one can excite a nondi verging whistler wave, 
i.e., a ducted whistler. The latter corresponds to a field 
aligned density hump with N ~ 10 -10 . Furthermore, the 
whistler solitons with density dips and humps can occur 
depending on the consideration of the frequency regime 
as well as the magnetic field strength and/ or the particle 
density. 

We ought to mention that our present investigation 
on the nonlinear propagation of EM whistlers might 
play an important role in studies of beat-wave parti- 
cle accelerators 43[ as well as in the problem of radio- 
frequency electron-cyclotron-resonance heating [44j of 
plasmas where the driver, instead of being a laser, is a 
whistler wave. 
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